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Model Independent Methods for Determining B (T(5S') 
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We describe a method that provides a model independent measurement of the Bs fraction in 
T(5S) resonance decays, fs, using the relative rates of like-sign versus opposite sign dileptons; the 
like-sign leptons result from B^ and Bs mixing. In addition, we show that determining the rates of 
single, double and triple Dg mesons provides an alternative way of finding fs- 

PACS numbers: 13.25.Hw, 13.66.Bc 



The T(5S) resonance has long been thought to be a 
source of Bs mesons, since it is massive enough to pro- 
duce BsBs pairs [1]. Recently, the CLEO collabora- 
tion established the presence of Bs mesons and made 
a model dependent measurement that the Bs fraction, 
fs, of T(5S') decays is (16.0±2.6±5.8)%, using a the- 
oretical estimate of (92±11}% for the inclusive branch- 
ing ratio B{Bs —>■ DsX) [3L ExpHcit Bs final states 
have also been reconstructed [3| . These results have been 
confirmed by Belle 4]. Precision measurements of Bs 
branching fractions are one of the most important goals 
of such studies. For example, it has been suggested that 
measurements of B{Bs -^ Dg Dg') lead to a deter- 
mination of the lifetime difference between CP-I- and CP- 
eigenstates [5j. To measure these, however, it is imper- 
ative that the number of Bs mesons produced in T{5S) 
decays be well known. The purpose of this paper is to de- 
scribe techniques that can be used to produce accurate, 
model independent measurements of fs- (Note that we 
are not considering T(5S) decays to modes that do not 
contain B mesons.) 

Both B and Bs mesons can result from T{5S) decays. 
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surement of Bs production [lO|. We must also account 
for B(i mixing and thus need to know the composition of 
the BB states, since the level of mixing will depend on 
the Charge Conjugation (C) state of the B''^B system. 

Hadrons containing 6-flavor are produced in pairs in 
current experiments. If the B and B are not in an eigen- 
statc of C, then the probability for a -6° to decay as a 
B , integrated over time, is given by 
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where the mixing parameter x = Am/T. 

The ratio R of mixed events to unmixed events is given 

by 
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N^o^o is equal to 
is equal to x^ + (1-a)^- 



The possible final states 
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For incoherent states then NgoB" 
2x(l - X) and N^„^o + NgO^„ is equal to x^ + (1 

For Bs mesons the mixing oscillations are so fast that 
R equals 1, regardless of the C parity state of the BsBs 
pair. For the B'^, however, R does depend on the initial 



configuration then [11 1 
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where the mixing parameter Xd = Amd/^ and is well 
measured as 0.775±0.008 [13]. 

For C even configurations we have 



B^*^B TT, BBtttt, and Bg Bs , since there is not suf- B^B state. If the ijiitial pair of B mesons is in a C odd 
ficient energy to produce an extra pion. 

Our first method for determining fs requires the mea- 
surement of like-sign versus opposite-sign dileptons. High 
momentum leptons from B (or Bs) decays reflect the 
flavor of the parent; positive leptons result from B de- 
cays, while negative leptons arise from B decays. We 
will assume here that the minimum lepton momentum 
requirement is large enough so that contamination from 
the decay sequence B -^ DX, D — » Yiv is negligible, or 
suitable corrections can be appfied [6|- 

This technique relies on the fact that ^5 mixing oscilla- 
tions are very rapid compared to Bd- The mass difference 
for B^ mesons is Arud = 0.509 ± 0.0005 ps~^, while for 
Bs mesons, Attic is limited at 90% confidence level to 
be >16.6 ps~^ [71, y, and recently measured by CDF to 
be 17.3lto;^8 ± O-O'' PS"^ @1- Bs mixing, when the me- 
son decays semileptonically, produces a relatively large 
number of like-sign dilepton events that allows a mea- 
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The states containing charged B pairs do not con- 
tribute to mixing. We expect that there are an equal 
number of such states as neutral non-strange i3's, since 
the mass difference is small and Coulomb corrections are 
small even at the T(4S) [H. The B^W and B*^'B'^ are 
both produced in an L equals one state; thus the B^B'^ 
pair is in a negative C parity state. The B^B'^^ state, 




FIG. 1: Ratio (R) of like-sign to opposite sign dileptons as a 
function of a; = Am/F for C odd (dashed), C even (dotted) 
and incoherent (sohd) B^B or BsBs states. 
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FIG. 2: Functions for like-sign dileptons (rising curves) and 
opposite-sign dileptons (falling curves) for C odd (dashed), C 
even (dotted) and incoherent (solid) B°B or BsBs states. 



(from B'^B^ or B^* B^) on the other hand, is in a posi- 
tive C parity state. 

Next we consider BBii and BBtttt states. The 
B'^B^tttt states can be incoherent, or coherent if there 
is fixed angular momentum between the pions and be- 
tween the pions and the B^s . R for an incoherent state 
is given in terms of x as (lj| 
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The variation of the R functions with x is shown in Fig.[T] 
Note that i?°_B~7r+ and B^B+tt" final states contain 
only one B^ and mix according to Eq. [l] leading to a 
mixing rate given by Eq. [31 and thus would be treated 
as an odd C state for our purpose. The other states can 
also be classified in terms of C eigenstates or incoherent 
states. 

CLEO has made the first measurement of B meson 
production at the T(5S) ^. JThey find_that B*"B*° 
production is largest with B*^B^+B^B*'^ being about 
1/3 its rate. Although limits on other decay channels are 
not small, we will assume that more data will allow the 
exact composition of the B decays at the T(5S) to be 
established. 

We now consider like-sign dilepton production coming 
from neutral S's. The yield from B^B pairs is given by 
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where iVss is the number of T(5S') events containing a 
BB pair above continuum background, and Ba-si is the 
B'^ semileptonic branching ratio. The Di^±±(x) func- 
tions for the neutral B pairs in C=± eigenstates, or being 
incoherent, are listed in Table HI and plotted in Fig. [51 
The fractions in each of these three states are denoted by 
fi; they include the BB'k{'k) final states [la|. We also list 



and plot the Di^±z^(x) functions for opposite sign dilep- 
tons. The fraction of neutral non-strange B mesons, fd, 
is assumed to be equal to that for charged B mesons [17 1 . 
Since the sum of neutral, charged, and strange B mesons 
is unity, we have 



fu + fd + fs = 1 

'^fd + fs = 1, and therefore 
fd = fu = (l-/s)/2. 

C state Di-±±{x) Di_±zp 
x^ {2 + x') 



(7) 



Odd 



Even 



Incoherent 



2(l-ha;2) 2(l + a;2) 

x'^{3 + x'^) {2 + x^ +x^) 
2(l-ha;2)2 2(H-x2)2 

x'^ + 2a^ 2 + 2^2 + a;4 

2(l+x2)2 2(l+x2)2 



TABLE I: Functions for like-sign and opposite-sign dileptons. 



A similar set of expressions exist for like-sign leptons 
from Bs decays 
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where Bs-si is the semileptonic branching ratio and 
D±±{xs) is the function that characterizes the dilepton 
rate and, in principle, depends on whether or not BsBs 
is in an even or odd eigenstate. The function form is 
identical to the D±±{x) functions listed in Table HI but 
incoherent states are not allowed since there isn't enough 
energy at the T(5S') to produce an additional pion. Sim- 
ilarly, the opposite-sign functions for Bs, D±^{xs), are 



identical in form with the D±z^{x) functions. Note that 
all these functions are normalized so that they go to a 
value of 0.5 as x (or xs) gets large, reflecting the fact 
that the mixing probability goes to its maximum value 
of 50%. 

The rate of opposite-sign dileptons is given by 



iV+- + A^-H 



N5S [fsBl_,,iD±^{x,) 



(9) 



l-fs 



^l.;E/»^-±t(^) 



1-/5 



b: 



-sl\ 



The last term is due to charged i?„ decays which do not 
mix. We form the ratio of like-sign to opposite-sign dilep- 
tons and divide through by the charged B semileptonic 
branching ratio B'^_gi. We replace the resulting ratios of 
semileptonic branching fractions by the lifetime ratios, 
Bd-si/Bu-si = t{B")/t{B+) = f+o, and Bssi/Busi = 
t{Bs)It{B^) = f+5. Since the CDF measurement of 
Ams when multiplied by the measured lifetime gives 
Xs «26.6 [9|, we substitute D±±{xs) = D±^{xs)—0.5. 
Denoting 
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the resulting equation is 

fsfls + (1 - fs)flo E^ /.A-±±(X) 

^ fsfls + (1 - fsno E. /.A-±t(^) + (1 - fs) 
Solving for fs gives 

Sd - P 



where 



Sd = fir 



^' {p - l)fls + Sd-p' 



^fiDi-±±{x) - p^fiDi^±^{x) 



(11) 
(12) 



(13) 



For practical application, we can use Eq. [7] to deter- 
mine fd and fu and then use well measured i?° and B^ 
branching ratios to normalize the Bs rates. 

We can estimate the error in fs by taking the C odd 
contribution as 75%, the C even component as 25%, no 
incoherent B^B contribution and fs equal to 16% from 
the CLEO model dependent determination that used 
B{Ds -^ 07r+) = (4.4 ± 0.5)%, which is an average be- 
tween the PDG value 7] and a recent BaBar measure- 
ment [l8|. In this case p=0.25. Taking into account the 
semileptonic branching ratio (10.5%), the fraction of high 
momentum leptons above the minimum Icpton momen- 
tum cut (1/3), and the lepton efficiency (0.8), we esti- 
mate that an error of ±4% on fs can be achieved with 
30 fb~^ of data. (We find that the fractional error in fs 
is about twice the fractional error in p.) Interestingly, a 
new preliminary value of B{Ds -^ (pn^) = (3.5 ± 0.4)% 



based on CLEO-c data [l9| raises fs to 21%, changes 
p to 0.29, and consequently improves the sensitivity to 
about ±3%. Thus, this method can lead to a model in- 
dependent determination of fs once the C states of the 
B'^ decays are experimentally determined more precisely. 

We next discuss another method for finding the Bs 
fraction that uses B mixing as input, but is not the main 
ingredient. Rather, we make use of the fact that Dg 
mesons are produced in almost all Bs decays, while they 
are produced only at the 10% level in B decays. We 
define B = B{B -^ DsXY. This rate has already been 
measured in T(4S) decays [2,ll|- One complication is that 
the Bs often decays into a DgDg pair, which is ignorable 
in B decays. To take into account the possibilities of 
the Bs decaying into one or two Ds mesons, we define 
Si = B{Bs -^ DsX), and S2 = B{Bs -^ D+D^X). 
Here 5*1 is the rate for the Bs to decay into one and only 
one Ds meson, and S2 is the rate for the Bs to decay into 
a Ds meson pair. The inclusive rate B{Bs -^ DsX) — 
Si + 25*2. (Note that the probability that a Bs does not 
decay into a Ds meson is given by 1 — Si — 6*2.) 

In this analysis there are two sources of Ds production, 
the Bs again with fraction fs, and B decays, where we 
take the mixture of i?° and B'^ events to have the same 
Ds yields as on the T(4S'). Since we do not expect the 
charge of the B to effect the Ds rates, this method is 
insensitive to their relative contribution. Thus, the B 
fraction is taken as (l-/s). 

Consider now the production of multiple Ds candi- 
dates in single T(5S') events, taking into account mix- 
ing. T(5S') decays can produce events with 4 Ds mesons, 
when both Bs mesons decay into DgDgX. We denote 
as TV"' ' the observed number of such 4 Ds events 
and note that this rate is not effected by mixing. N^^^ 
refers to events with 3 observed Ds, whose charge sum 
is either +1 or -1, and whose rate is also not effected by 
mixing. N'^ denotes finding an event with oppositely 
charge DgDg mesons, whose rate is changed by both 
Bs and B mixing, so we introduce the parameters /^j^ 
and fmix, where f^i^= 1/2 and /„.„ equals average mix- 
ing rate over the C even, odd, incoherent mixtures and 
charged B decays defined above. N^^ denotes a DgDg 
or DgDg pair, while N^ indicates the detection of a 
single Dg meson. Here, the single rate is inclusive of all 
double, triple and quadruple rates, etc.. The resulting 
equations relating the observed numbers to the branch- 
ing ratios and fs are 



N+~+-/e^N5s 

N^^^/e^N,s 

N^^/e^N5s 
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Slfs 

{2SiS2 + ASl)fs 



(14) 
(15) 

[(1 - fL.)Sl + 2(1 - fL.)SiS2 (16) 
+2S2{l-Si-S2)+^Sl]fs 

[fLM + 2SiS2) + 2Si)]fs (17) 



+ (1 - fs)B fmix 

N^/eN^s = 2(5i + 2S2)fs + 2(1 - fs)B , (18) 

where e indicates the detection efficiency of a single Ds ; 
it is the sum of the branching ratio times efficiency for 
each decay mode that is used. We ignore two smah cor- 
rections; first of aU we take the rate for a single B me- 
son to produce a DgDg pair to be neghgibly small, and 
secondly we don't account properly for the charge cor- 
relations resulting from the small production of "wrong- 
sign" Ds production from Bs decays that can occur via 
a b —> u transitions. The latter effects only equations 
containing /^„ ^. 

We expect that e could be made as large as 10% by 
including many modes. We use 5*1 = 0.8 and S2 = 0.1, 
which allows estimates of the various rates. We do not 
expect to be able to observe a significant quadruple Ds 
rate. On, the other hand, we estimate that in a 50 fb^^ 
data sample, there would be ~500 observable triple Ds 
events allowing the use of Eq. [151 Thus, we would have 
four equations, in this case Egs. llStfTSl relating our three 
unknowns: /s. Si and 5*2. The best values for the un- 
knowns can be obtained by using a constrained fit to find 
the solution. Note that Eq. [T6]can be simplified as 



N^'F/e^N; 



-2S2 



5S - [(1 - fmix)Sl - 2/; 
2Sl]fs + {\-fs){l-Lm 



,SiS2 (19) 

)B2 . 



Should not enough triple Ds events be found in the 
data sample, Eqs. fTMTSl could be solved for the three 
unknowns. Another possibility, if precise measurements 
on the B^ C parity contributions are not available, is to 
add Eqs. [TBI and [T71 since the resulting equation is not a 
function of either B^ or Bs mixing. In this case, however, 
measurement of the triple Ds rate is necessary. 

In conclusion, knowledge of the number of Bs mesons 
is essential for all branching fraction determinations at 
the T(5S'). We present a model independent method of 
determining the fraction oi Bs mesons at the T(55'), and 

hence B ( T(5S') — » Bg Bs ) , using the complete mixing 
of the Bs via dileptons. The amount of luminosity re- 
quired to make an accurate measurement of the Bs frac- 
tion will depend on the actual compositions of the B^ 
final states, but it is likely to require several tens of fb"^. 
We also suggest another technique using the counting of 
single, double and triple Ds mesons in T{5S) events 21 1. 
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In addition, to take precisely into account any "wrong- 
sign" B decays, e. g. B~ -^ DfX, we advocate mea- 
suring fmix on the T(4S') using the ratio of like-sign to 
opposite-sign Ds events. 

It may also be possible to determine B and B~ produc- 
tion cross-sections by reconstructing exclusive decays. All 
possible modes B^*'B , B^*^ B vr, and BBtv-k must be 
included. The sum of these cross-sections with the Bs 
yield found using the techniques described above should 



equal the total T(5S) cross-section. Any short fall would containing b and b quarks, 

be indicative of decays of the T(5S) to final states not 



